Stabilizer Quantum Error Correction with Qubus Computation 
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In this paper we investigate stabilizer quantum error correction codes using controlled phase 
rotations of strong coherent probe states. We explicitly describe two methods to measure the Pauli 
operators which generate the stabilizer group of a quantum code. First, we show how to measure 
a Pauli operator acting on physical qubits using a single coherent state with large average photon 
number, displacement operations, and photon detection. Second, we show how to measure the 
stabilizer operators fault-tolerantly by the deterministic preparation of coherent cat states along 
with one-bit teleportations between a qubit-like encoding of coherent states and physical qubits. 

PACS numbers: 03.67.Pp, 03.67.Lx, 42.65.-k 



The question of which physical system is best suited for 
quantum information processing is still open, each imple- 
mentation proposal having strengths and weaknesses. In 
some systems (such as optics) it is difficult to make qubits 
interact, so that the two-qubit gates needed for univer- 
sal computation are difficult to implement. One scheme, 
proposed by Gottesman et. al. [H, circumvents the need 
to make qubits interact directly by using a modified tele- 
portation protocol. A generalization of this leads to the 
cluster state proposal of Raussendorf et. al. where 
a large entangled state is prepared offline, and computa- 
tion is performed by a sequence of single qubit measure- 
ments which depend on the outcomes of previous mea- 
surements. A different scheme that bypasses the need for 
qubits to interact directly was proposed by Nemoto et. 
al. 0, HI- This scheme shows that, by inducing a phase 
on a large coherent state bus mode which depends on the 
logical state of the physical qubits, one can implement a 
near deterministic CNOT gate between the qubits. Co- 
herent states are particularly useful because of the ease 
with which they may be produced, e.g. with lasers or 
Bose-Einstein condensates. Further developments have 
shown more direct methods to perform two-qubit gates 
with bus modes, termed qubus computation [5j. 

If qubus computation is to be seriously considered for 
physical implementation, a full analysis of the propaga- 
tion of errors should be undertaken. The starting point 
for these considerations is whether we can perform quan- 
tum error correction (QEC) on qubits efficiently. In par- 
ticular, can one measure the syndromes for a given stabi- 
lizer code directly with controlled rotations (CRs) and 
strong coherent probe beams? Recent work by Yam- 
aguchi et. al. @ demonstrates how to measure the syn- 
dromes for some stabilizer codes using these tools. They 
show that the stabilizers for the three bit-flip code can be 
measured directly with CRs and a single strong coherent 
bus mode. The stabilizers for Shor's 9-qubit code can 
also be measured, showing that it is possible to correct 
for any error on a single qubit in an encoded block. 



The purpose of this paper is to generalise the results of 
Yamaguchi et. al and demonstrate how CRs can be used 
to implement quantum error correction with any possible 
stabilizer code. We will describe two schemes to measure 
the syndromes of an arbitrary weight n Pauli operator, 
using the stabilizer operators of the seven qubit code as 
a concrete example for each one of these schemes. The 
first scheme uses a single strong coherent probe beam, 
a quadratic number of CRs, a linear number of coherent 
displacements, and a photon number measurement. This 
scheme can be modified to use homodyne measurement 
at the cost of a slightly larger number of CRs and coher- 
ent displacements. The second scheme we describe is a 
fault-tolerant approach to the measurement of the Pauli 
operators, which requires a linear number of strong co- 
herent pulses, CRs and detectors. Although we focus on 
the 7 qubit code - which has stabilizer generators with 
weight 4 - for each of these schemes we describe how to 
generalise to Pauli operators of weight n. 

Background — In Q it was shown that the stabilizers 
for the 3-qubit bit-flip code (|0)->|000), UM 111 )) c °uld 
be measured with the parity gate depicted in Fig. [1^,. 
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FIG. 1: (a) Circuit to measure the parity of two qubits. The 
CRs are ±0. (b) Two parity gates combined to measure the Pauli 
operators ZZI and IZZ. 

It can be seen that this circuit is a parity gate when 
we consider its effect on the input state j^in) = (ao|00) + 
ai|01) + a 2 |10) +a 3 |ll))|a), where |a) is a coherent bus 
mode. The effect of the CRs is to apply a phase to the 
coherent beam if our data qubit is |1) and leave it alone 
otherwise: (a|0) +6|l))|a) -»■ a\0)\a) +b\l)\ae w ). Before 
the detector D 2 in Fig. QJi, the state \tpi n ) is (a |00) + 
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a 3 |ll))|a) +ai|01)|ae- ie ) + a 2 |10)|ae ie ). When D 2 is a 
homodync detection along the x-quadrature we are able 
to distinguish a |00) +a 3 |ll) from ai|01) +a 2 |10), since 
a homodyne measurement of \a) along the x-quadrature 
is equivalent to the projection (x\a). That is, \ae ±l6 ) are 
indistinguishable when we homodyne detect along the x- 
quadrature. This is the basis of the CNOT shown in 

With two parity gates we can measure the Pauli oper- 
ators ZZI and IZZ. That is, one parity gate is applied 
to qubits 1 and 2 to measure ZZI while the second par- 
ity gate is applied to qubits 2 and 3 to measure IZZ, 
as shown in Fig. [TJd. The state before the application 
of the parity gates is \ipi n ) = (co|000) + ci|lll))|a)|a). 
There are four cases to consider: no error, |^>; n ); an error 
on qubit 1, XII\ifj in ); an error on qubit 2, IXI\ip ia }] an 
error on qubit 3, IIX\ip in ). We can see what the effect 
of a bit flip error on each of the modes is by considering 
the state |a6c)|a)|a), where a, b, c S {0, 1}. Directly be- 
fore homodyne detection in Fig. [Tb |a6c)|a)|a) becomes 
\abc)\ae l ^ a ~ h ^ B )\ae l ^ b ^ 6 ) . When we measure the probe 
states to be \ae ±lm8 ) \ae ±m6 ), where m,n € {0,±1}, we 
know whether there was no error (m, n = 0) or a one 
bit flip error, the location of the bit flip also being iden- 
tified by the values of m and n. Similar methods can 
be applied to measure the stabilizer operators for Shor's 
9-qubit code. The natural question that arises is: can 
we use techniques similar to those above to measure the 
syndromes for an arbitrary stabilizer code? 

larger codes — As a concrete example, consider the 
[[7, 1,3]] stabilizer code 0]. This code can correct a sin- 
gle arbitrary quantum error in any of the 7 qubits, and 
it has been used extensively in studies of fault-tolerance 
in quantum computers due to the fact that it allows for 
simple constructions of fault-tolerant encoded gates [1]. 
In order to detect which error has corrupted the data, 
one must measure six multiqubit Pauli operators which, 
up to qubit permutations and local unitaries, are equiv- 
alent to the Pauli operator ZZZZ, or the measurement 
of only the parity of 4 qubits. For an arbitrary stabilizer 
code, various multiqubit Pauli operator must be mea- 
sured, each of which is always equivalent to a measure- 
ment of only the parity of a subset of qubits, thus it 
is sufficient to consider only multiqubit parity measure- 
ments in order to perform quantum error correction with 
stabilizer codes. 

Single Coherent State Pulse — In order to measure 
ZZZZ with CRs, we can start with the encoded state 
(co|0l) + ci|li))|a), and design a circuit that gives us 
\a±) when there was no error (even parity) and \a 2 ) when 
there was an error (odd parity), where a± =/= a 2 - 

Ideally we would want to do this with just one coherent 
probe beam, four CRs and a single homodyne detection, 
following a direct analogy with the circuit depicted in 
Fig. [TJa. However this is not possible. The best we can 
do is have some even states go to \a) and the rest go to 
\ae ±2lB ) while the odd states go to \ae ±lB ). The circuit 



that performs this is shown in Fig. [2^i. 




FIG. 2: (a) First attempt at using CRs to measure the parity of 
four qubits with a single strong coherent probe, (b) Position of 
phase space peaks corresponding to the state of the coherent probe 
beam. Yellow circles correspond to odd states while red and black 
circles correspond to even states. 

Notice that in phase space we would have five points 
- three for the even states (\a), \ae ±2lB }) and two for the 
odd states (\ae ±l0 )) - as can be seen in Fig. [2b- If we were 
to homodyne detect the probe beam at this stage, we 
would partially decode our encoded state co|0l) +ci|1_l) 
since we can distinguish the state \a) from \ae ±2q,e ). The 
problem now becomes determining what operations must 
be done before we homodyne detect so that we only dis- 
tinguish between states of different parity in the first four 
qubits, and nothing more. It turns out that either homo- 
dyne or photon number detection can be used, depending 
on the operations applied before the measurement. 

Photon number measurement — If we incorporate dis- 
placements along with Fig.[2K we can take the five points 
in phase space to just three. Displacements of a state 
can be easily implemented by mixing the state with a 
large coherent state on a weak beam splitter, the size 
of the coherent states amplitude and beam splitter re- 
flectivity deciding the displacement [H . If we have three 
displacements and three applications of Fig. [5^, as in 
Fig.Ek,, we find that |odd) -> | - 4asin 2 (0/2)(2 sin 2 (0) + 
cos(0))) and |evcn) -> | ± 2asin 2 (0)(2cos(0) - 1)), as 
depicted in Fig. [3b- The displacements that accom- 
plish this are D{fi x ) = D{-4a cos 2 (6/2) (2 cos(0) - l)), 
D(J3 2 ) = D(a(l + 2cos(0) + 2cos(30))) and D(J3 3 ) = 
£>(a(cos(20) - cos(30) - cos(0) - 1). 
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FIG. 3: (a) CRs used to measure the parity of four qubits with 
a photon number detection, (b) Position of phase space peaks 
corresponding to the state of the coherent probe beam. 

Notice that the red and black circle in the Fig. [5b are 
equidistant from the p-axis. We can thus perform a pho- 
ton number measurement on the probe beam to deter- 
mine whether we had an odd or even state. In order for 
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a photon number measurement to distinguish the odd 
from even states we require a.9 2 3> l/y/3. 

We can use this method to measure the parity for a 
state of any size. If we have n qubits then we can have 
at best n+ 1 points in phase space, using the 9, —9, 9, —9 
pattern for the CRs shown in Fig. [2^. Using displace- 
ments and a photon number detector we are able to mea- 
sure the parity. In general, if n is even we need n — 1 
displacements and n 2 — n CRs with a photon number 
measurement. When n is odd, after the application of 
the circuit analogous to Fig. [2K of size n, we will have 
the point \ae l ( 2n ~ 1 * >e ) in phase space without the point 
\ae~ lty2n ~ 1 ^ B ) . So we need an extra displacement to move 
the non-symmetric point. If n is odd we need n displace- 
ments and n 2 CRs with a photon number measurement. 

For this method to work we need the use of a num- 
ber discriminating photo-detector. In practice it is well 
known that homodyne detection is much more precise 
than number discriminating photo-detectors. For this 
reason, we describe how to measure a Pauli operator us- 
ing homodyne detection. 

Homodyne detection - - Consider the ZZZZ case 
again. After applying Fig. [2^ we have five points in phase 
space. Ideally we want |a) and \ae ±2t8 ) to become one 
point in phase space, say i?i+«i?2, and |ae ±J0 } to become 
one point, say R3 + iR±. If this was possible then homo- 
dyne detection could be used. This can be done with five 
displacements and six applications of Fig. [2^, requiring 
10 simultaneous equations to be solved. Without loss of 
generality we set R2 = R4 = Ri — 0. The equations to 
be solved are e 4i/3fJ 'A + e 3i/3< 'B + e 2lfiB C + e m D + E = 



equivalent to measuring the parity of 4 qubits and noth- 
ing else), we would use the circuit shown in Fig.[4ji. 



-i/38 



dp 



e 5ip8^ w h ere ao ±2 = and a±i = R3 



After solving these equations we find that A, B, C, 
D and E scale as -R z /9 i . We are free to choose the 
distance between the origin and R3 to be arbitrarily large, 
at the expense of using arbitrarily large displacements. 
We can also use the above method to distinguish the 
parity of any given state of n qubits. If we have n qubits 
we have n+1 points in phase space, using a circuit similar 
to Fig. [2^,. In order to distinguish the parity we need n+1 
displacements and n(n + 2) CRs. 

Fault- Tolerance — These two methods to measure 
weight n Pauli operators cannot be used for fault-tolerant 
quantum computation. If there is an error on the coher- 
ent probe mode during one of the CRs, say photon loss, 
it would be transferred to a phase error in each of the 
physical qubits it interacts with afterwards - that is, a 
single fault can cause a number of errors which is greater 
than the number of errors the code can correct. For this 
reason we now look at measuring the syndromes of sta- 
bilizers fault-tolerantly. 

Shor [ljj first described how to fault-tolerantly mea- 
sure the generators of the stabilizer group of a quantum 
error correcting code using ancilla GHZ states HO}®" + 
|l}® n )/V2, CNOT's and Hadamards. For example, in 
order to measure the Pauli operator ZZZZ (which is 




FIG. 4: (a) Circuit for the measurement of the parity of four 
qubits Q3 . (b) Same circuit modified to use coherent states and 
CRs. 



To fault-tolerantly measure the stabilizer group gen- 
erators of a QEC with CRs we make three modifica- 
tions to Fig. 0k.. First, instead of using |0) and |1) for 
the ancilla, we use the coherent states |a) and \ae l6 ), 
respectively. In that case, the ancilla GHZ state be- 
comes (\a)® n + \ae i9 )® n )/V2. Second, we replace the 
CNOT's with CRs, which will cause a phase shift if 
the physical state is |1) and do nothing otherwise, i.e. 
|l)|a) -> \l)\ae~' W ) and \l)\ae t6 ) -» |l)|a). We also 
need to replace the Hadamards with some quantum oper- 
ation H which will perform the mapping H\a) m (|ct) + 
\ae i9 ))/V2 and H\ae w ) « (|a) - \ae l9 ))/V2. Third, we 
replace the qubit measurements with some sort of optical 
measurement that distinguishes between \a) and |ae ±l9 ) 
but not between \ae l6 ) and \ae~ ie ) - this is what we call 
Z measurement. This new circuit is depicted in Fig.[4)D. 

The Z measurements can be performed directly by 
homodyne detection, or by displacements followed by 
photon counting detectors - in both cases, using tech- 
niques outlined earlier in this paper. What remains to 
be specified is the preparation of the coherent cat state 
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ton . 



\ae l6 )® n ) / 'y/2 and the implementation of the H 



operation. One solution for the cat state preparation is 
to use one bit teleportations ll| which translate states 
from the |0)/|1) basis to the \ot) /\ae ±l6 ) basis. Prepa- 
ration of the cat state is done by using the one bit tele- 
portation in Fig. to prepare \\fna) + \y/nae lB ) from 
the state (JO) + |l))/-\/2 and the coherent state \^/na), 
and then sending this state into an n-port symmetric 
beam-splitter [12|, |l3|. In principle, we are required 



to correct the state before the beam-splitter by apply- 
ing the transformation Z such that Z\a) ~ \a) while 
Z\ae l8 ) « — \ae l6 ). However, we can avoid explicitly ap- 
plying this transformation by keeping track of this neces- 
sary correction - what is called the Pauli frame T|- and 
compensating for it in subsequent measurements. Simi- 
larly, to perform the H (the approximate Hadamard on 
coherent state logic) we first use Fig. [SJd to telcport the 
quantum state from the bus to a qubit, then perform the 
Hadamard transformation and finally teleport back to 
the coherent state logic using the circuit shown in Fig. [5^,. 
These teleportations, when performed back-to-back to 
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FIG. 5: Approximate one bit teleportation protocols llll using 
CRs. The measurement dependent Pauli corrections need not be 
performed, as discussed in the body of the paper. 



teleport a qubit state to another qubit, can also be used 
as leakage reduction units to reduce leakage faults to reg- 
ular faults [l5j |. 

The resources required to measure a weight n Pauli 
operator are 3n + 1 CRs, n + 1 ancillary qubits, In Z 
measurements and n + 1 qubit measurements. 

Noisy ancillas — If the probability of error at each 
gate is bounded by e, transversal operations and encod- 
ing can ensure that the probability of an uncorrectable 
error is 0(e 2 ) instead of 0(e). An error during cat state 
preparation may lead to correlated X-like errors in the 
cat state with probability O(e), which can lead to uncor- 
rectable errors in the encoded data during the measure- 
ment of the Pauli operator, thus defeating the purpose 
of encoding the data for fault-tolerant quantum compu- 
tation. In order to avoid this, one can verify the in- 
tegrity of the cat state via non-destructive state mea- 
surement 



16l . Il7j ] . When using CRs and coherent beam 



probes, this translates to preparing an extra copy of the 
cat state, which remains in coherent state logic, inter- 
acting with the qubit GHZ state transversally with con- 
trolled —6 rotations, and Z measuring each mode of the 
ancillary cat state. By performing classical error cor- 
rection on the measurement outcomes, one can deduce 
the locations of A-like errors in cither the GHZ state or 
the ancillary cat state. If the data is encoded in a code 
that can correct a single error, repeating this procedure 
with another ancillary cat state allows for the inference 
of which locations in the qubit GHZ state have X errors 
with high enough probability to ensure uncorrectable er- 
rors are only introduced into the data with probability 
0(e 2 ) [la ], so that Pauli measurements with a verified 
ancilla can be used for fault-tolerant quantum computa- 
tion. Overall, the overhead for each attempt of measuring 
a weight n Pauli operator consists of 2(n + 1) CRs, 2 an- 
cillary qubit preparations and measurements, and 2n Z 
measurements. 

Z-likc errors (including dephasing of coherent super- 
positions, one of the consequences of photon loss in the 
CRs) do not lead to errors in the encoded data, just er- 
rors in the outcome of the Pauli operator measurement. 
If error correction is to be performed, the Pauli operator 
measurement must be repeated 3 times, and a majority 
vote of the outcomes is taken, in order to ensure that the 
measurement outcome is reliable [l6| . 

Some of the systematic errors in the probe beams, such 
as phase rotation or attenuation (also consequences of 



photon loss in the CRs), can be partially compensated 
for by additional linear-optics elements and by adjusting 
the Z measurements individually to minimize additional 
X errors. Moreover, errors in the transversal operations 
during the preparation of the cat state are independent, 
and thus do not need special consideration during this 
verification stage - they do contribute to e, however, and 
are thus crucial for fault-tolerance threshold calculations. 

Discussion — We have shown two schemes to mea- 
sure the syndromes of an arbitrary weight n Pauli opera- 
tor. The first scheme uses a single strong coherent probe 
beam, a quadratic number of CRs, a linear number of co- 
herent displacements, and a photon number or homodync 
measurement - however, this scheme is not fault-tolerant. 
The second scheme we described is fault-tolerant, and the 
amount of resources scales linearly with the weight of the 
Pauli operator. This demonstrates how it is in principle 
possible to perform general fault-tolerant quantum com- 
putation in the qubus architecture. It is worth noting 
that we could have easily used controlled displacements 
in the place of CRs in the methods presented here. 
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